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THE NASTRAN HYDROELASTIC ANALYZER 


James B. Mason 
Test and Evaluation Division 


SUMMARY 

An overview of the NASTRAN Hydroelastic Analyzer com- 
puter program is given. Included is a theoretical derivation show- 
ing general assumptions and procedures used in obtaining the finite 
element model of the hydroelastic problem. Program implementa- 
tion of the theory is discussed, and an Appendix showing a compar- 
ison between NASTF '*N and theoretical results is attached. 
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THE NASTRAN HYDROELASTIC ANALYZER 


♦ The motions of liquids and gases in various containers has been of interest for 
many years,* One particular area where these effects are of great importance 
is in the analysis and design of missile and space vehicle structures. In these 

* type structures, where lightweight-thin- walled construction is typical, the con- 
tainers cannot be considered rigid and coupling of the various fluid motions with 
the elastic deformations of the tank must be considered. This so-called hydro- 
elastic problem and, in particular, its solution using the NASTRAN computer 
program is the subject of this report. 

The NASTRAN (NASA STRUCTURAL ANALYSIS) computer program is a general 
purpose finite element program developed for the analysis of large and/or com- 
plex structures. It was decided that the addition of certain fluid capability to this 
program would result in a useful hydroelastic analyzer for the solution of many 
important problems. Other advantages arising from this approach would be that 
many automated features already existing in NASTRAN would apply to the resulting 
hydroelastic analyzer. For example, since the existing structural program had 
the capability of treating control systems, the resulting hydroelastic analyzer 
could include the coupled interactions of the fluid, the structure, and the control 
system for problems where these influences are of concern. 

The addition of the hydroelastic capability was accomplished by the MacNeal- 
Schwendler Corporation under contract to the Goddard Space Flight Center. This 
capability is presently operational on the IBM 360-95 computer. 

It is the purpose of this report to present the general capabilities, limitations, 
and characteristics of this new hydroelastic analyzer. To accomplish this, a 
theoretical derivation is given first which presents general assumptions and 
procedures used in obtaining a finite element statement of the hydroelastic prob- 
lem. The derivation is not limited by topology considerations and discretization 
is symbolically indicated with no associated finite element property specification. 

The program implementation of the theory is then discussed. Topological limita- 
, tions and fluid element properties employed in the program are given. General 

input/output data and user selected solution requests are reviewed. The report 
concludes with a summary of the capabilities and limitations associated with the 
- hydroelastic analyzer. 

* For example, see Abramson, H. N., "The Dynamic Behavior of Liquids in Moving Containers," 

NASA SP- 106. 
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An Appendix is included showing a comparison of NASTRAN and theoretical 
results for several sample problems. Gases and liquids with free surfaces are 
treated here, and eigenvalue and forced response solutions are obtained. 

The basis for the analytical developments reported herein is the work by the 

MaeNeal Schwendler Corporation and, in particular, the derivations by Dr. R. t 

H. MaeNeal and D. N. Herting, see references (1) and (2). The theoretical 

developments in this report, however, present an alternative formulation of 

many aspects of the hydroelastic problem. It is hoped that this reformulation v 

will provide the prospective user with both an overflew of the program and a 

document which can serve as an aid when used in conjunction with the above 

references. 


THEORETICAL DEVELOPMENT 

The problem to be treated is that of a compressible fluid in an elastic container. 
Only small motions from the undistrubed state of the coupled fluid- structure 
system are considered. We desire the pressure tq be the basic unknown of the 
fluid and the deformation to be the basic unknown for the structure. 

In this section, the governing differential equations describing the pressure dis- 
tribution in the fluid continium are derived. A variational statement equivalent 
to these governing equations and boundary conditions is then obtained. From 
this variational statement, the matrix equations of the finite element fluid model 
are indicated. The final set of finite element matrix equations describing the 
hydroelastic problem are then obtained by combining the equations describing 
the fluid and structural systems. 


Fluid Description - Differential Equations 


We assume that the mass density p d of the disturbed fluid never varies greatly 
from its undisturbed value P . Thus, we may write 


P d = p (1 + s) (1) 

where s = s (x, y, z, t) is the condensation and is small compared to unity. 

We also assume that disturbances propagate through the fluid by an adiabatic 
process. With this assumption, we relate the pressure p to the density by the 
state equation 
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p = p (p d ) 


( 2 ) 


Expanding (2) about the undisturbed condition and using (X), we obtain 


■ p = p + JLEL* s p + . . higher order t 


'° T ^I 


erms 


( 3 ) 


where p 0 is the undisturbed pressure. For an adiabatic process 


P<) y 

P=— P\ 

P v 


( 4 ) 


where y is the ratio of specific heats at constant pressure to constant volume. 
Substituting (4) into (3) and using (1) yields 


p =p Q + — 7 s p = p Q + a 2* s p = p 0 + a*-<p d - p) 

P y 


( 5 ) 


= constant + a* p d 


In the above, the speed of sound a 0 through the undisturbed fluid has been used 
since 


y — p\ rX) $ y^'Cl + (y - 1) s] Syf? = a 2 (6) 

p y P p 0 

Neglecting nonlinear and viscous terms in the equations of motion and assuming 
that the fluid velocities are small compared with a 0 , Euler's equations take the 
form* 

* Note that ILL? = ( • ) 

9 t V ' 


3 



(?) 


where u is a small displacement of the fluid from a fixed reference point and 
only gravitational body forces p d g are assumed acting. Substituting (1) into (7) 
and neglecting products of small terms, results in the equations of motion 


p 


• • 
u 


-Vp-pg 


•* 


( 8 ) 


The continuity equation can be expressed as 


P A + u = 0 


(9) 


Using (5), equation (9) ma^ be written as 



a 


2 

o 


p + V* p d u - 0 


( 10 ) 


Substituting (1) into (10) and neglecting products of small terms results in the 
equation of continuity in the form 


-L p + p V* u -- 0 (H) 

a o 


*■■4 

The particle velocity u may be eliminated from (8) and (11) by taking the first 
time derivative of (11) and the divergence of (8), then the three dimensional 
wave equation 


P a; 


v-ivp * v-g 

p 


( 12 ) 
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in terms of pressure is obtained. The quantity p 0 a* is equal to the adiabatic 
bulk modulus B of the fluid. Thus, equation (12) can be written in the form 



- V* 1 Vp s v-g 

p 


(13) 


lh the following developments the effects of the gravity gradient are neglected 
and the wave equation takes the final form 


i p -v-ivp = 0 

B p 


(14) 


The boundary conditions to be considered are of two types: 


p= prescribed on boundary surfaces S 2 


(15a) 


i 2-? - - q_ on boundary surface S, 
p o n M * 


(15b) 


Here d p / 3 n is the derivative of p in the direction of the outer unit normal n to 
the boundary surface and q n is a generalized normal flux per unit surface area. 
It will be shown subsequently that condition (15a) corresponds to fluid free 
surface points in the absence of a gravity field while (15b) corresponds to the 
conditions at both a free surface when gravity is present and points on a fluid- 
structure boundary. 


Fluid Description - Finite Element 


While the variational statement of a problem is not required for the application 
of the finite element method, this approach does provide a clarity and consistency 
to the finite element formulation and is adopted here. 
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From the results of variational calculus, see references (3) and (4), it is known 
that for the functional* 


H (p) 



(volume) 



q n P d S 
(surface) 


(16) 

r 


to attain a minimum, then the pressure must satisfy the Euler equation (14) 
through the volume and the conditions (15) on the boundary. Thus, equation (16) 
can be viewed as an alternative statement of the fluid problem. Taking the 
variation of (16) and considering p as an invariant, we have 




6 p ■+• b (JL Vp -^p\ld V 


I 


so that 


A P & p d V + JiVp-SVpdV=~J 


q H p dS = 0 ( 17 > 


q n opdS 


(18) 


represents the variational statement of the problem for the continuous fluid 
system. We note here the similarity of (18) with the principle of virtual work 
for dynamic loading in a continuous structural system. From this viewpoint, 
the first term on the left of (18) can be associated with the virtual work of 
"gene ralized inertia forces" (l/B)p , the second term with the "generalized 
strain energy" of the fluid volume, and the right hand side of the equation with 
the virtual work of "distributed generalized surface forces" q n . 

The discretization of the fluid system in terms of a finite number of degrees of 
freedom p ■ is obtained from the variational statement (18) using standard 
procedures (see references (3) and (5)). To acomplish this, we relate the 
pressure p in the continuous system to the degrees of freedom p. - { p t p • • • } 
by the matrix equation 


* Note that 


p is considered invariant and p is net constrained on S r 
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(19) 


P = Ca] (p.) 

where [a] is a function of position. Similarly, the pressure gradient may be re- 
lated to (pj) by the matrix expression* 

vp = {Vp} = [b] (p.) (20) 

From (19) and (20), we obtain the virtual quantities 

S p = [a] (Sp.) ( 21a ) 

S Vp = 8 {Vp} = [b] { a Pj ) (21b) 

Substituting (19) through (21) into (18) and noting that the virtual pressures are 
arbitrary, results in the matrix equations of motion 

[M F 3 (Pp + [K F Hp.} = {I' > (22) 

for the idealized fluid system. In (22) we have 


on = f [a'] T i [a] d V (23) 

Jv 


* Note that the form of [a] and [b] is determined by the character of the spatial pressure 
distribution assumed for the finite fluid elements. Here we symbolically represent these 
matrices and defer discussion of the fluid element properties until the next section. 




( 24 ) 



where the superscript F refers to the fluid. Examination of (23) and (24) shows 
that the "generalized fluid mass matrix" [M F ] includes fluid compressibility 
effects while fluid mass effects are included in the "generalized fluid stiffness 
matrix" Ck f ]. 

The vector {I } represents "concentrated generalized surface forces" due to the 
distributed generalized surface forces q n . We observe that the forces I. do 
work through their conjugate generalized degrees of freedom p . and that I . 
is nonzero only on those degrees of freedom associated with the S x boundary 
where q n is defined. 

Two general types of boundary conditions are considered in the hydroelastic 
problem; they are: 

(1) Fluid-structure boundary 

(2) Fluid free surface boundary 

where the second condition is associated with gravity wave effects of liquids in 
containers. With these considerations, it is convenient to reorder and partition 
(22) into the form 



where the subscripts in the above are fluid degree of freedom identifications 
defined as 
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i = degrees of freedom associated with the fluid interior 
f = degrees of freedom associated with the fluid free surface boundary 
s = degrees of freedom associated with the fluid-structure boundary 

Expressions for the concentrated generalized surface forces I, and If applied, 
respectively, to the fluid-structure and free surface degrees of freedom will 
now be obtained. Examination of (25) shows that this requires knowledge of the 
generalized surface force (or flux) at these boundaries. 

At a fluid-structure boundary, uo fluid may penetrate the structural boundary 
and no gaps may occur between the fluid and the structure. For the nonviscous 
fluid considered, this means that the fluid acceleration component normal to 
the boundary must equal the acceleration of the structural boundary normal to 
itself. From (8) and the acceleration compatibility condition, we may write 


1 Bp 

p Tn 


•*« -4 A 

u* - g • n 
n w 


(27) 


at the boundary. In the above: 


l 

P 


B P _ 


i 

1 Vp-n 

P 


n = outward unit normal at the fluid-structure boundary 

u*s m normal component of the structural acceleration at the fluid- structure 
boundary 

Comparing (27) with (15b) we see that the generalized normal flux at the fluid- 
structure boundary is 


# 


q 


n 


ir 


n 


-* A 

+ g *n 


( 28 ) 


Substituting (28) into (25) and neglecting the static gravity effects (since all 
motions are measured from static equilibrium) results in the expression 
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(29) 



[a,]T U SdS 


for the concentrated forces applied to the fluid-structure freedoms p . S is 
the fluid- structure boundary area portion of S* . * 

The boundary condition problem at a free surface is rather unique. This arises 
from the fact that while the boundary condition at the free surface is known, i.e. 


p = pressure at the free surface -0 


(30) 


since measurement is made from static equilibrium, the location of this boundary 
is not. The method adopted for the program stems from the Eulerian descrip- 
tion used in fluid mechanics where the analysis views a point fixed in space, 
see sketch below 



To determine the pressure p at a point in space originally 1 °ated on the free 
surface at static equilibrium, we expand in a Taylor’s series about p on the 
displaced surface. Thus, we obtain 


P = p- u*Vp + 


- u-Vp 


( 31 ) 


since p = 0 and higher order terms have been neglected. Using (8) in (31) and 
neglecting higher order terms, we obtain 
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(32) 


p = p u *(u + g) $ p g-u 


Assuming that g is positive in the direction of the outward normal to the original 
free surface permits (32) to be written in the form 


P = P 8 uj 


(33) 


where 

g s magnitude of the gravity vector 

uF 5 fluid free surface displacement component parallel to the gravity vector. 
Differentiating (33) twice with respect to time yields 


»• 


P 


= pgu 


F 

n 


(34) 


By a procedure like that leading to (28) we may obtain the generalized normal 
flux at the free surface as 


q n = U J + g-n 


(35) 


Substituting (34) into (35) yields 


q = JL + 5* 

n P g 


(36) 


Substituting (36) into (25) and neglecting static effects results in the expression 




-I 


[aJTjL dS 


P g 


(37) 
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for the concentrated forces applied to the free surface freedoms p f . S f is the 
free surface boundary area portion of S x . 

In the absence of a gravity field we see from (32) that the free surface boundary 
condition becomes 


p ~ 0 (38) 


Comparing (38) with (15a) and the variational statement of the fluid problem, we 
see that (38) is an "essential boundary condition" while those represented by 
(29) and (37) are "natural boundary conditions." Cases represented by (38), 
where the pressure is prescribed to be zero over an area S 2 , are treated auto- 
matically in the program by the application of constraints requiring 



(39) 


The results of (29) and (37) are now to be put into a form suitable for use in a 
general matrix statement of the hydroelastic problem. To accomplish this, the 
character of the finite fluid elements defining [a s 3 and f 3 must be known and 
assumptions and numerical procedures for the integrations defined. In addition, 
for [I s ] the independent degrees of freedom and the type of finite elements used 
for representation of the idealized structure must be considered* Complete 
details of the procedure are presented in references (1) and (2). 

For the present, it will be satisfactory to assume that the discretization process 
associated with (29) and (37) has been accomplished. Symbolically indicating 
this for the concentrated forces applied to the fluid degrees of freedom p at the 
fluid- structure boundary, we write (29) in the form 


«.}=- [A1 T (u*s J 


(40) 
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where 


ujj s normal structural displacement components of the finite element 
structural model at structural degrees of freedom on the fluid- 
structural boundary associated with fluid degrees of freedom p g *. 

[A] T = rectangular matrix describing the interconnection of the fluid and 
structure. This matrix reflects the structural degrees of freedom 
used in the analysis and the structural boundary areas associated 
with these freedoms. 

Similarly, the concentrated forces applied to the fluid degrees of freedom p f at 
the free surface can be indicated by writing (37) as 


n f > = cm] (p f > 


(41) 


where 


[m]= square matrix describing the free surface effects. This matrix in- 
cludes gravity and fluid density effects. 

Equations (40) and (41) may now be included in the general matrix statement of 
the fluid problem. Noting that (41) has the form of a generalized fluid mass 
coefficient, we may write (26) in the form 


Mfi | Mf f ; m?; 


t 

\* a 


1 m 

fa— » 

L* 

< i 


- - 


i * 

j T ~ ~ 

Mff j «f f ♦ “ | Mf. 
“f. ! Mf. \ m;. 

< 

• e . 

Pf 

*• . 

P. 

► + 

i 1 

Kfr ; Kf f ; K f. 

Kf, j KJ 4 I K*. 

— | 1 — 




— * • - 


r *> 

Pi 


0 

<pf 

► = < 

0 

p« 

w - 


- [A ] 1 {Uy 

w * 


(42) 


* Note that the subscripts s and s refer to fluid and structural degrees of freedom, respectively, 
on the fluid structure boundary. 'There is no requirement that s =*s and, in fact, this will usually 
not be the case. 
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Structure Description - Finite Element 

The finite element displacement formulation of the structural dynamics problem 
is well known and is not repeated here, see references (2), (3), and (5) for de- 
tails . For our purposes, we express the Lagrangian matrix statement of the 
discretized structural mechanics problem in the partitioned form 



where 

u® - structural displacement components normal to the fluid-structure 
boundary. 

= structural displacement components tangent to the fluid-structure 
r boundary 

u r = structural displacement components for the remainder of the struc- 
ture, i.e., those degrees of freedom describing the structure which 

are not included in [u s ! u s 3 . 

n r i V 

[ = structural mass matrix 

[K s ] = structural stiffness matrix 

P, ; , P t _, P r = applied external loads (measured from static equilibrium) which do 
work through the conjugate structural degrees of freedom u®_ , 
us , us ; respectively. s 

t- r 

p n-’ p u “ l° ac * s Produced by the fluid (measured from static equilibrium) 

which do work through the conjugate structural degrees of freedom 
u5 » u t » > respectively. 

5 S 

Equation (43) differs from that of the standard structural mechanics finite ele- 
ment formulation only in the additional loading terms P F _ and P F U due to the 
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coupling at the fluid-structure boundary. Since viscous effects in the fluid are 
neglected*, we have the structural boundary conditions 


{P F }=0 

V 



(44a) 

(44b) 


where S y is the area associated with structural degree of freedom u® r and p 
is the fluid pressure adjacent to the displaced fluid-structural boundary. Differ- 
ences in the Eulerian description of the fluid and the Lagrangran description of 
the structure requires a treatment similar to that leading to (32) in order to 
define p , see sketch below. 



Expanding about p and neglecting higher order terms, we have 


P = P- /?g*u s (45) 


Writing g and the structural displacement vector u s in terms of their components 
normal and tangential to the original surface and substituting into (44b), we 
obtain 


* See references (1) and (2) for suggested methods of including these effects. 
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(46) 


• pF \ 

n r‘ 



[p - f, g t uf 



d S 


Equation (46) must now be put into a form suitable for use in a general matrix 
statement of the hydroelastic problem. As was done in presenting (40) and (41), 
we shall again symbolically indicate the discretization process associated with 
(46) and write this in the form 


KJ = CA] {p,} - [K n? 




(47) 


where [A3 has been introduced before and [K] is a rectangular matrix reflecting 
the differences in Fulerian and Lagrangian description. The [K] matrix is a 
function of gravity. 

Equation (47) may now be included in the general matrix statement of the struc- 
tural problem. Noting that the [K] matrix in (47) has the form of a structural 
stiffness coefficient, we may write (43) in the form 






r._ ^ 

MS 

m® . ! 

M s 


u s 

n s " t ; 

"T ff J 

n, r 


ii 

M s , ' 

Mf 1 

Mf 

<4 

iif 

"I * «' * 
1 


r 



i 

-i 

- — 


— 

M s ! 

M*, r 

M rr _ 


ii 

_ n » r 

T 



- r - 


K s - K 

; K s , - K„. 

i 

K s 1 


kl 

n ? n . 

; n T'„ 
i 

i 

i 

n » r 


", 

K s „ 

J KS t 

“r 

i 

KS 

4 

u s 

n r l . 

| 

i 

t, r 




-4 — 

i 

— 



K s 

! K?. 

i 

K s 


U 

L "» r 

* t; r 

i 

rr _ 


r 

V. J 


P n, * <P,> 


r < 48 > 



Coupled Fluid- Structure Description - Finite Element 

By combining the matrix equations (42) and (48) representing, respectively, the 
fluid and the structure; we obtain the final system of equations describing the 
hydroelastic problem. 
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i-:4 r ; -fw 
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k ' **• j t"* j 


Equation (49) is symbolic of the general unsymmetric system of finite element 
equations describing the hydroelastic problem. A discussion of the problem 
types considered by the NASTRAN hydroelastic analyzer and the various possible 
problem solution types is given in the following section. 


GENERAL DESCRIPTION OF NASTRAN HYDROELASTIC ANALYZER 

The previous derivations have been rather general in nature and applicable to 
any fluid and/or structural topology consistent with the limiting assumptions. 

The intent has been to provide a general overview of the theory employed by the 
program for the solution of the hydroelastic problem. While general methods 
have been indicated, the details of implementation are not given since this would 
be beyond the scope of the present report. The details of this implementation 
are fully discussed in references (1) and (2). 

It is the purpose of the present section to briefly specify the problem type 
treated by NASTRAN, indicate the general input/output data associated with 
the program, and describe the fluid finite element properties used. 

The program has been written to treat the case of axisymmetric fluid volumes 
(axisymmetric in both geometry and fluid properties) held in elastic containers. 
However, neither the elastic properties of the structure or the motions of the 
coupled fluid- structure system need be axisymmetric. Several sample fluid 
topologies which could be treated by the program are shown in Figure 1. 

The NASTRAN formulation of the hydroelastic problem uses the Fourier coef- 
ficients of the pressure in an axisymmetric coordinate system as the independent 
degrees of freedom for the fluid, and uses the structural grid point displacements 
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Figure I Examples of Permissible Fluid Topologies 


and rotations as the independent degrees of freedom for the structure. All 
NASTRAN structural elements with the exception of the axisymmetric elements 
(i.e. conic shell, doubly curved shell, and solid of revolution) can be used for 
the hydroelastic analysis. 

Input data to the hydroelastic analyzer associated with the Quid* can be briefly 
summarized as consisting of the following information: 

1. Cross-sectional coordinates of fluid points (fluid circles) — RINGFL cards 

2. Connection information and physical properties (, and B) for Quid ele- 
ments connecting fluid points - CFLUIDi cards 

3. Identification of structural grid points - GRIDB cards 




Structural input data is basically unchanged from that of the standard NASTRAN as it exists 
for structural purposes. 
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4. Sequential list of fluid {joints on free surfaces and flu id- structure 
boundaries - FSLIST and BDYLIST cards 

5. The harmonic conqjonents of fluid pressure to be considered - AXIF 
cards 

6. The magnitude of the gravity vector - AXIF cards 

7. Identification of any radial planes of symmetry - FLSYM cards 

8. The location of {joints (r, z, p) on the free surface at which normal dis- 
placements are desired for output - FREFPT cards 

9. The location of {joints (r, z, 4) in the fluid at which values of the pres- 
sure are desired for output - PRESPT cards 

10. Direct input matrix terms coupling fluid degrees of freedom with euih 
other or with any structural degree of freedom or extra points - DMIAX 


The above input data will be discussed in more detail below. The axisym metric 
fluid volume is described by the analyst using a number of concentric circular 
rings (equivently ’’fluid points”) distributed through the fluid volume and con- 
nected by fluid elements having the properties of bulk modules and density, see 
Figure 2. The user defines the location of the various rings in the fluid volume 
by specifying the location of the intersection point (fluid point) cf the ring with 
an arbitrary transverse r-z plane containing the axis of symmetry. The fluid 
point locations are given in either a cylindrical or spherical coordinate system 
on the RINGFL input data cards while CFLFDi input cards are used to specify 
fluid element connection data. 
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cards. 


f 


Three types of fluid elements are provided: 

1. CFLUID3 - Triangular element defining an annular volume of triangular 
cross- section. This element connects three RINGFL points. 


z A 


RINGFL POINT (TYPICAL) 




2, CFLUID4 - Quadrilateral element defining an annular volume of quadri- 
lateral cross-section. This element connects four RINGFL points. 



3. CFLUID2 — Center element defining a fluid disk. This element is re- 
quired to connect a fluid volume to the axis of fluid symmetry. The 
center element connects two RINGFL points and defines a fluid volume 
bounded by two parallel planes perpendicular to the fluid axis of sym- 
metry and a conical outer boundary. 
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Figure 2 shows a sample fluid model utilizing the various fluid elements. 

The fluid element properties defining the character of the finite element matrices 
(i.e. of [a] and [b] of the assembled system) are obtained by expanding the pres- 
sure in a Fourier series with respect to the circumferential coordinate 4 > , thus 


N N 

p (r, z, <p) = P n cos n <j> + ^T" 1 P n * sin n 0 
n* 0 n* 1 


(50) 


where P n and P n * are functions of the axial (z) and radial (r) coordinates. As- 
suming, in turn, a functional relation over the r,z cross-section of the fluid ele 
ment in terms of generalized coordinates q? and q " * , we can write 


p n _ p n ( r , z) = [L n (r, z)j {qj} n >0 
pn* = pn* (Ti z ) = jjn z )] - ( q n* } n ;> 1 


(51) 


Here the matrix CL n j reflects the assumed pressure distribution over the ele- 
ment cross-section, and the generalized coordinates are chosen so that the 
values of P n and P n * match the values of the pressure coefficients at the fluid 
points located at the corners of the element. Thus, if PJ is the value of P n at 
the i th fluid grid point, we may write 


{qy> = toy • <py} 

(qp = [HJ p ] {Pp 


(52) 


Substituting (52) and (51) into (50) results in an expression for p in terms of the 
fluid degrees of freedom PJ and PJ* . This relation is equivalent, on an ele- 
ment level, to equation (19) with p . designating the fluid element freedoms. The 
element properties, i.e. generalized mass and stiffness, can then be obtained by 
the procedure indicated by equations (23) and (24). 
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From the above discussion, we see that each fluid point has a number of har- 
monic pressure degrees of freedom associated with it. In this respect we note 
the similarity with the six degrees of freedom usually associated with a struc- 
tural grid point in NASTRAN. 

The triangular fluid element properties are derived under the assumption that 
the pressure varies linearly with position over the cross-section described by 
the triangle. For this case 


CL n (r, z)] = [1 • r i zl- 


and 


(53) 



I 

1 I r 2 
I 



(54) 


where r i9 (i = 1, 2, 3) are the coordinates of the three connecting fluid points. 

The quadrilateral fluid element properties are obtained by averaging four tri- 
angular elements; see sketch below. 



The properties of the center element are derived using an assumed pressure 
distribution obtained from an asymptotic solution of the wave equation as 
r - 0. For this case 

[L n (r, z)] = [(r) n j (r) n z] (55) 
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[H n ] = 
1 qp 


(r,) n 


(O n 


(r,) n z , 


( r 2 ) n z 2 


(56) 


where r. , z. (i = 1, 2) are the coordinates of the two connecting fluid points. 

The effects of spatial variations in fluid properties, e.g. the case of oil on water, 
can be included in the analysis since fluid properties are specified at the ele- 
ment level. 

In addition to the RINGFL and CFLUIDi cards already described, additional in- 
put Bulk Data Cards provided for the hydroelastic problems are briefly reviewed: 

1. The AXIF card describes the fluid coordinate system, the gravity 
parameter parallel to the axis of fluid symmetry, and a user chosen list 
of fluid harmonics to be used for the particular problem. 

2. The FSLIST and BDYLIST cards describe the fluid boundary. The 
FSLIST card lists the fluid points which lie on a free surface while the 
BDYLIST card lists fluid points on a fluid-structure boundary. Fluid 
boundary points not included on either of the above lists are assumed 
to lie on a rigid boundary. 

3. The PRESPT and FREEPT data cards specify special fluid output 
points. The PRESFI card specifies locations for pressure output (i.e., 
in addition to harmonic pressure output). The FREEPT card is used to 
define locations on the free surface boundary for surface displacement 
output (see equation (33)). 

4. The FLSYM data card allows the user to optionally model portions of the 
problem with planes of symmetry (symmetric and/or antisymmetric) 
containing the fluid axis of symmetry. 

5. The GRIDB card identifies structural grid points for hydroelastic prob- 
lems. This card identifies the specified structural grid point with a 
particular RINGFL fluid point. (It should be noted that the r, z coordi- 
nates of each structural grid point on the fluid-structure boundary must 
coincide with the r, z coordinates of one of the fluid points on the inter- 
face.) 
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6. The DMIAX card allows the user to input special purpose matrices for 
his particular problem. For example, surface friction effects can be 
included by methods introduced in references (1) and (2). 

From the given input data, the program automatically assembles the hydroelas- 
tic problem, including the effects of free surfaces and fluid-structure boundaries, 
represented by (49). Various solutions associated with (49) may then be re- 
quested from the twelve Rigid Formats available within NASTRAN. 

The available Rigid Formats are: 

#1 - Static Analysis 
#2 - Static Analysis with Inertia Relief 
#3 - Normal Mode Analysis 
#4 - Static Analysis with Differential Stiffness 
#5 - Buckling Analysis 
#6 - Piecewise Linear Analysis 
#7 - Direct Complex Eigenvalue Analysis 
#8 - Direct Frequency and Random Response 
#9 - Direct Transient Response 
#10 - Modal Complex Eigenvalue Anr lysis 
#11 - Modal Frequency and Random Response 
#12 - Modal Transient Response 

While the fluid effects have been programmed so as to be compatible with all of 
these, the following characteristics of the hydroelastic problem cause restric- 
tions on the type of non- trivial solutions which may be obtained: 

1. The formulation is for dynamic perturbations from static equilibrium. 

2. The system of equations (49) describing the general hydroelastic problem 
is unsymmetric. Since the first six Rigid Formats are restricted to the 
use of symmetric matrices, the program automatically ignores the fluid- 
structure boundary. Thus, for these Rigid Formats the program solves 
the problem of a fluid in a rigid container with an optional free surface 
(equivalent to (42) with (u®-} = {0} )and an uncoupled elastic structure 
with no fluid present (equivalent to (43) with {P£_} = { P*^ = {0}). 

3. The only direct means of exciting a fluid is through the fluid-structure 
boundary since no means is provided for the direct input of applied loads 
on the fluid. 

With the above restrictions in mind, the suggested Rigid Formats are outlined 
below: 
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Rigid Format #3 — The uncoupled modes of a fluid in a rigid container and 
a structure with no fluid may be found. 

Rigid Format #7 - The coupled modes of the fluid and structure may be 
found. 

Rigid Formats #8 and #9 - External loads may be applied to the structural 
grid points only. For the Direct Transient Analysis, initial conditions may 
only be given for the structural degrees of freedom since the fluid is as- 
sumed to be initially at static equilibrium. 

Rigid Formats #10, #11, and #12 - The practicality of these is limited since 
the modal coordinates used to formulate the dynamic matrices are the nor- 
mal modes of the fluid and structure solved as uncoupled systems. For 
widely spaced fluid and structure frequency ranges, however, some utility 
may exist for these formats. 

Upon solution of the particular problem, the program outputs both fluid and 
structural data. The available structural data is unchanged from the present 
NASTRAN structural program. Fluid output data which may be directly re- 
quested includes fluid harmonic pressures, fluid pressures (PRESPT), and free 
surface displacements (FREEPT). First and second time derivatives of fluid 
pressure and free surface displacement may also be obtained. Printed values 
for the fluid may include real and complex values. Both x-y Plot and random 
analyses capability are available for FREEPT and PRESPT points. 


CONCLUDING REMARKS 

A derivation leading to a general finite element statement of the hydroelastic 
problem has been presented and a description of the NASTRAN hydroelastic 
analyzer based on this derivation has been given. Some of the more important 
characteristics of the hydroelastic analyzer are summarized below. 

Limiting restrictions: 

1. Only small dynamic perturbations from static equilibrium are considered. 

2. The fluid must have axisymmetric geometry and properties. 

3. Fluid disturbances propagate by an adiabatic process. 


4. The fluid is nonviscous. 



5. Gravity gradient effects are neglected. 

6. I1‘ gravity is present, the gravity vector must be parallel to the axis of 
fluid symmetry and any free surface must lie in a plane perpendicular 
to this axis. 

7 . The application of external loads is allowed only on structural grid 
points. 

8. Axisymmetric structural elements cannot be used. 

9. The effects of variable fluid density and gravity field are not completely 
compatible. If the effects of either gravity or variable density are small, 
they may be used together but second order errors may result. 


Capabilities and Characteristics: 

1. Motions of both coupled and uncoupled fluid- structural systems may be 
treated. 

2. Rigid and/or elastic structural boundaries can be included. 

3. Multiple boundary effects may be analyzed. 

4. Effects of variable fluid density and compressibilities may be included. 

5. Fluid free- surface effects can be treated. 

6. Degrees of freedom of the fluid are the harmonic pressure coefficients. 

7 . Degrees of freedom of the structure are the grid point displacements and 
rotations. 

8. Fluid behavior is measured from a fixed reference point. 

9. Symmetric and antisymmetric capability exists to take advantage of 
symmetry planes. 

10. Fluid description is generally compatible with all NASTRAN system 
capability. 
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An experienced user of the program using available NASTRAN features (e.g. 
Direct Matrix Inputs, Alters, Multi- Point- Constraints, etc.) can circumvent 
many of the restrictions in the above summary lists. Several examples of this 
are briefly discussed here. 

As has been mentioned previously, viscous effects can be included by the input 
of a damping matrix to (49) for the inclusion of viscous wall friction. Addi- 
tionally, sources and/or sinks in fluids without free surfaces could be treated 
by the addition of ’’generalized body forces” at the source/sink points. Also, 
one and two dimensional Laplacian fluid elements, obtained from existing 
NASTRAN scalar and membrane elements by procedures explained in reference 
(6), could be used in conjunction with the CFLUIDi elements for treating certain 
fluid problems*. 

The NASTRAN hydroelastic analyzer should provide a powerful analytical tool. 
The present program capability of control system modeling can be included to 
yield a description of coupled fluid- structure- control system problems. The 
direct application of the program to many acoustics problems is obvious. 


* It should be noted that one and two-dimensional Laplacian elements could be used "stand-alone" 
to treat many problems if care is taken in the choice of generalized mass and boundary terms. 
Also, with this approach, velocity potential, steam functions, pressure, etc. could be taken as 
the fluid freedoms. • 
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APPENDIX A 


This material is taken directly from the NASTRAN 
Demonstration Problem Manual and shows a com- 
parison between hydroelastic analyzer and theoretical 
results for several sample problems. 



RIGID F0RMAT No, 3, Real Eigenvalue Analysis 
Vibration of a Compressible Gas in a Rigid Spherical Tank (3*2) 


A. D escription 

This problem demonstrates a compressible gas in a rigid spherical container. In NASTRAN a 
rigid boundary is the default for the fluid and, as such, no elements or boundary lists are 
necessary to model the container. 

Aside from the NASTRAN bulk data cards currently implemented, this problem demonstrates the 
use of the hydroelastic data cards: AXIF, CFLUID2, CFLUID3, and RINGFL. 

The lowest mode 'frequencies and their mode shapes for n ■ 0, 1 and 2 are analyzed where n Is 
the Fourier harmonic number. Only the cosine series Is analyzed. 

B. Model 

1 . Parameters 

R - 10.0 m (Radius of sphere) 

p ■ 1.0 x 10" 2 3 Kg/m 3 (Mass density of fluid) 

B * 1.0 x 10 3 Newton/m 3 (Bulk modulus of fluid) 

2. Figure 1 and 2 show the finite element model. The last 3 digits of the RINGFL identifica- 
tion number correspond approximately to the angle, 0, from the polar axis along a 
meridian. 


C. Theory 

From Reference 18, the pressure in the cylinder Is proportional to a scries of functions: 


4 n,m 


Jmxi (X) „ 

— i pj* (cos e) cos n<j> , 

ST m 


n < m 

m » 0,1 ,2 


0 ) 


where: 


Q n m Pressure coefficient for each mode 

•I fin 

(i) k 

X Nondimenslonal radius * r 

a 

&>mk Natural frequency for the ktb mode number and mth radial number in radians 
per second 

0^ Bessel function of the first kind 
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p 


radius 


a speed of sound In the gas 

pJJ associated Legendre functions 

6 merldlnal angle 

$ circumferential angle 

n harmonic number 

m number of radial node lines 

The solution for X and hence Is found by the use of the boundary condition that the flow 
through the container Is zero. 



where R Is the outer radius . 

This results in zero frequency for the first root. Multiple roots for other modes can be 
seen in Table 1. The finite element model assumes different pressure distributions In the two 
angular directions which causes the difference in frequencies. 

D. Results 

Table 1 and Figure 3 summarize the NASTRAN and analytic results for the lowest nonzero root 
in each harmonic. Table 1 lists the theoretical natural frequencies, the NASTRAN frequencies, the 
percent error In frequency, and the maximum percent error In pressure at the wall as compared to 
the maximum value. Figure 3 shows the distribution of the harmonic pressure at the wall versus 
the meridinal angle. The theoretical pressure distributions correspond to the Legendre functions 
P 0 {cos e), (cos 9), and Pj (cos 6) which are proportional to cos e, sin e, and sin e respec- 
tively. 
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Table 1. Comparison of NASTRAN and analytical results 



Natural Frequency (Hertz) 

Pressure 

Analytical 

NASTRAN 

% Error 

Rvini 

0 

33.1279 

33.2383 

0,33 

1.19 

1 

33.1279 

33.2060 

0.24 

0.47 

2 

53.1915 

53.3352 

0.27 

0.91 
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Figure 1. Gas filled rigid spherical tank model. 


CFLUID2 and CFLUID3 
Identification Numbers 
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Figure 2. Gas filled rigid spherical tank model. 
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Figure 3. Pressure at tank wall - first finite modes. 





RIGID FORMAT No. 3, Real Eigenvalue Analysis 
Vibration of a Liquid In a Half Filled Rigid Sphere (3-3) 


A. Oescri ptlon 

The model Is similar to Demonstration Problem No. 3-2 except that a hemispherical fluid model 
with a free surface Is analyzed. Additional cards demonstrated are the free surface list (FSLIST) 
and free surface points (FREEPT), The effective gravity for the fluid Is found on the AXIF card. 
The fluid Is considered Incompressible. 

The lowest three eigenvalues and eigenvectors for the cosine and sine series of n » 1 are 
analyzed, where n Is the harmonic order. 


B. Input 


1. Parameters 


9 

R 


10.0 ft/sec‘ 
10.0 ft 


p • 1.255014 lb-sec 2 * /ft 4 


B 


(Gravity) 

(Radius of hemisphere) 

(Fluid mass density) 

(Bulk modulus of fluid. Incompressible) 


2. Figure 1 shows the finite element model. 

C. Results 

Reference 17 gives the derivations and analytical results. In particular, the parameters 
used In the reference are: 


e » 0 (half-filled sphere) , 


A ■ ~~ (dimensionless eigenvalue) . 

9 


0 ) 


Table 2 of Reference 17 lists the eigenvalues, A-j, A 2 » and A 3 for the first three modes. 
Figure 13 of Reference 17 shows the mode shapes. 

The analytic and NASTRAN results are compared in Table 1. The frequencies are listed and 
the resulting percentage errors are given. The maximum percent error of the surface displacement, 
relative to the largest displacement. Is tabulated for each mode. 
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Tht fret surface displacements nay be obtained by th* aquation 


where p la the praaaura at the free surface recorded In the HAST RAN output. Note that* tinea an 
Eultrlan reference frame Is used* the pressure at the original (undisturbed) surface Is equal to 
the gravity head produced by motions of the surface. Special FREEPT data cards could also have 
been used for output. Since the results are scaled for normalization anyway, the harmonic pres- 
sures may be used dlrectfy as displacements. 

Figure 2 Is a graph of the shape of the free surface for each distinct root. Both analytic 
and NASTRAN results are scaled to unit maximum displacements. Because the cosine series and the 
sine series produce Identical eigenvalues* the resulting eigenvectors may be linear combinations 
of both series. In other words the points of maximum displacement will not necessarily occur at 
$ ■ 0* or $ ■ 90°. Since the results are scaled* however, and the results at $ ■ 0 are propor- 
tional to the results at any other angle, the results at <p ■ 0 were used. 


Table 1. Comparison of natural frequencies and free surface mode shapes from the reference and 

NASTRAN. 


Mode 

Number 

Natural 

Frequency (Hertz) 

Mode Shape 


NASTRAN 

NASTRAN 
% Error 

Maximum 
% Error ,e 

1 

0.1991 

0.1983 

-0.1 

c < 1 % 


0.3678 

0.3691 

0.3 

e < 2.6% 


0.4684 

0.4766 

1.8 

e < 4 % 
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1090 2090 30 




Free Surface 



11 \ 10 ^ 



19 I TfU 


caui03 


CFLUID2 


419 

20175 20171 20167 


1 - 419 CaUIDi Identification Number 

1090 - 20175 RINGFL Identification Number 
(Where last 3 digits corresoond 
to the meridian angle and y 
leading digits correspond to / 
radial location) A- 


Rigid Container 


Figure 1. Rigid sphere half filled with a liquid. 
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Figure 2. Free surface mode shapes. 
3.3-4 (9/1/70) 


A- 10 


RIGID FORMAT Mo. 7, Complex Eigenvalue Analysis - Direct Formulation 
Complex Eigenvalue Analysis of a Gas-Filled Thin Elastic Cylinder (7-2) 


A. Description 

This problem demonstrates the use of symmetry to analyze specific harmonics of a fluid-filled 
structure. The problem to be solved consists of a cylindrical section filled with a compressible 
fluid. The end conditions for the cylinder and the fluid are two planes of antisymmetry* perpen- 
dicular to the axis. These end conditions correspond to the conditions that exist at periodic 
Intervals along a long* fluid-filled pipe vibrating In one of Its vibration modes. The antisym- 
metric boundary for the structure Is defined by constraining the motions which lie In the plane. 

An antisymmetric boundary for the fluid corresponds to zero pressure. This may be modeled* In 
NASTRAN. by defining the plane of antisymmetry as a free surface with zero gravity. 

The lowest natural frequencies and mode shapes for the third and fifth harmonics are analyzed 
separately. For the third harmonic* the structure Is defined by a section of a cylinder having an 
arc of 30 degrees or 1/12 of a circle. The fifth harmonic analysis uses a section having an arc 
of 18 degrees or 1/20 of a circle. The longitudinal edges* which were cut* are planes of symmetry 
and antisymmetry In order to model a quarter cosine wave length. 

The bulk data cards used are; AXIF* BDYLIST, CFLUID2, CFLUID4, C0RD2C, CQUAD1 * EIGC, FLSYM, 
FSLIST * GRIDB, MAT1 * PQUAD1 * RINGa, and SPC1. 

B. Input 

The finite element model for the third harmonic Is shown In Figures 1 and 2. Parameters used 
are: 

B 

p f 

p s 

E 

G 
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2.88 x 10 3 lb/1n 2 

(Bulk modulus of fluid) 

1.8 x 10' 2 lb-sec 2 /1n 4 

(Fluid mass density) 

6.0 x 10* 2 lb-sec 2 /1n 4 

(Structure mass density) 

1.6 x 10 5 lb/1n 2 

(Young's modulus for structure) 

6.0 x 10 4 lb/1n 2 

(Shear modulus for structure) 


A-ll 


10*0 Inch 
10.0 Inch 
0.01 Inch 


(Radius of cylinder) 
(Length of cylinder) 
(Thickness of cylinder) 


a ■ 
t ■ 
h ■ 

The model for the fifth harmonic It similar to the third harmonic model except that the angle 
covered by the structure Is 18* Instead t. 30*. This Is accomplished by simply removing the strife-, 
tural elements and boundary GRID8 points corresponding to the two right-hand layers of structure 
(between 18* and 30*). The FLSYM, FSLIST and 'PCI cards are changed to reflect the 1/20 symmetry. 

C. Theory 

The derivations and results for this problem are described In Reference 16. The results for 
various dimensionless parameters are listed. The particular parameters for the problem at hand 
are: 





P fl a 

0 ■ -gp ■ 0.0 ■ 


where n Is the ratio of fluid mass to structure mass, C Is the ratio of the wave velocity In the 
structure material to the wave velocity In the fluid. 0 Is the factor describing static pressuri- 
zation, P 0 . 

The basic assumptions for this analysis are: 

1. Thin shell theory Is used for the structure. The bending moment terms In the force 
equilibrium equations are Ignored In the results. 

2. The fluid is nonvlscous, Irrotatlonal , and small motions are onl; ro.^Id' d. 

This particular problem becomes relatively easy to solve since the moaa for the fluid 

In a rigid container and the modes of the structure with no enclosed fluid have the same spatial 
function at the interface. Each mode of the fluid Is excUed by only one mode of the structure 
and each mode of the structure Is excited by one mode of the fluid. The pressure In the fluid Is 
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assumed to be a series of functions 


P • P n e 1wt cosn* sin *5- Q n (r,w) , (1) 

where Q n it * lessel Function or a modified Vessel Function of the first kind. 

The characteristic shapes of the structure are a series of the form; 

u ■ A e* wt cos n$ sin ^ , (2) 

where u Is the displacement normal to the surface. The fundamental momentum equation for the 
fluid flow at the boundary Is: 

V(P(r))*e r • - p f u , (3) 

where e r Is a unit vector In the radial direction. 

The forces on the structure at the boundary are: 

3 2 F 

P(*i) » y - p s hu , (4) 

a az* 8 

where the function Fj Is defined by the differential equation on the surface: 

• t fi • < 5 > 

oZ 

The solution for Is obtained by assuming that 

Fj » 0 e 1wt cos n$ sln^r . (6) 

Combining Equations 1 through 6 results In the relationships: 
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Equation (7) Is a statement of the continuity of displacement. Equation (8) states the 
balance of the pressures. The above equations may be solved by Iterating on o. Reference 16 
provides solutions for m over a wide range of parameters. 

The analytic and NASTRAN eigenvalues are listed In Table 1. The corresponding errors In the 
eigenvalues are tabulated and the maximum errors In displacement at the container wall are given 
as the percentage of the maximum value. The cor..a1ner displacements in the radial direction at 
<p • 0.0 are compared In Figure 3. 


Table 1. Comparison of analytical and NASTRAN results. 


Harmonic 

Natural Frequency (Hertz) 

Mode Shape 

Analytical 

NASTRAN 

% Error 

Max. % Error 
In Radial OIsdI. 

3 

1.579 

1.595 

1.0 

» 

cT 

• 

o 

5 

1.011 

1 .049 

3.4 

0.5 
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Free Surface (Antisyniaetry) Anti syirr* try 

CFLUIO Identification Number - Fluid CQUA31 Identification Number 

Structure 





























































































































































Radial Displacement 
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Figure 3. Radial displacement for harmonic 3 and 5 normal inodes. 



RIGID FORMAT No. 9, Transient Analysis - Direct Formulation 
Transient Analysis of a Fluid-Filled Elastic Cylinder (9-3) 


A. Description 

The fluid-filled shell, used for analysis of the third harmonic. In Demonstration Problem 
No* 7-2 Is subjected to a step change In external pressure at t ■ 0 of the form 

p * p 0 sin y cos n$ 

The fluid Is assumed Incompressible In order to obtain an analytical solution with reasonable 
effort. The harmonic used Is n * 3. 

In addition to the cards of Demonstration Problem No. 7-2, DAREA, PRESPT, TL0AD2, and TSTEP 
cards are also used. Selected displacements and pressures are plotted against time. 


3. Input 

The finite element model Is shown in Figures 1 and 2. Parameters used are: 


B * 

09 

(Bulk modulus of fluid - Incompressible) 

p f ’ 

1 .# a W Z lb~sec 2 / 1 n 4 

(Fluid mass density) 

P s - 

7 2 4 

€.0 * >{/ lb-sec /In 

(Structure mass density) 

E « 

1.6 xlO 5 lb/fn 2 

(Young's modulus for structure) 

G - 

6.0 x 10 4 lb/ 1 n 2 

(Shear modulus for structure) 

a * 

10,0 inch 

(Radius of cylinder) 

% * 

10.0 Inch 

(Length of cylinder) 

n * 

0.01 Inch 

(Thickness of cylinder wall) 

p * 

K o 

2.0 

(Pressure load coefficient) 
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C. Theory 


The theory was derived with the aid of Reference 16 as In Demonstration Problem No. 7-2. 
Since the fluid Is Incompressible, It acts on the structure like a pure mass. Neglecting the 
bending stiffness, the equation of force on the structure Is: 


Ps 


« (m + ny ) w 



0 ) 


where: 

p # Is the leading pressure on the structure (positive outward), 
m ■ p f h Is the mass per area of the structure, 
irif Is the apparent mz* of lent fluid, 
w Is ihc ns'.mai displacement (positive outward) 


The function F Is defined by the equation. 



Eh 

a 


a 2 w 

’<■ A 


( 2 ) 


The spatial functions of pressure, displacement, and function F may be written in the form: 


P s * P 0 sin -j- cos n* 

772 

w ■ w Q sin-—- cos n4 
F ■ F Q s1n y cos n$ 

where p Q , w Q , and F Q are variables with respect to time only. 
Substituting Equations 3 Into Equation 2 we obtain: 

F . . Id (if Us 

0 

Substituting Equations 3 and 4 Into Equation 1 we obtain: 
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(3) 


(4) 
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The Incompressible fluid Is described by the differential equation: 

V 2 p • 0 * (6) 

Applying the appropriate boundary conditions to Equation 6 results In the pressure distribu- 
tion: 


p - p f Sin 32. cos(n$) I n (2f) . (7) 

where I n is the modified Bessel function of the first kind and p r Is an undetermined variable. The 
balance of pressure and flow at the boundary of the fluid, with no structural effects. Is described 
by the equations: 


"o * ■ p r *n (if) • < 8 > 

* -&L • (9) 

• r-*a 

Substituting Equations 3 and 7 into Equation 9 results In: 


->f 8 « * -f 1 ; (f) p . 


( 10 ) 


Eliminating with Equations 8 and 10 gives the expression for apparent mass, tn^: 


>.(¥) 


0f W O 


l n\ i f 


m f w o 


( 11 ) 


Substituting the expression for m^ from Equation 11 into Equation 5 results in a simple single 
degree of freedom system. When the applied leading pressure Is a step function at t * 9, 
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where 


and 


and 



* 



D. Resul ts 

A transient analysis was performed for the case n ■ 3 on the model and various displacements 
and pressures were output versus time up to one second. The theoretical frequency Is calculated to 
be 1.580 Hertz and the period Is 0.633 seconds. The displacements at two points on the structure 
(Point 91 Is located at $ - 0, z * 5.0; Point 94 Is located at <J> » 18°, z * 5.0) are plotted 
versus time In FI gore 3. 

The maximum error for the first full cycle occurs at the end of .<«e cycle. The ratio of the 
error to maximum displacement Is 4.7556. Changes In the time step used In the transient integration 
algorithm did not affect the accuracy to any great extent. The most probable causes for error were 
the mesh size of the model and the method used to apply the distributed load. The applied load was 
calculated by multiplying the pressure value at the point by an associated area. The "consistent 
method" of assuming a cubic polynomial displacement and Integrating would eliminate the extraneous 
response of higher modes. The tnethod chosen in this problem, however, Is typical of actual 
applications. 
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